One-dimensional itinerant ferromagnets with Heisenberg symmetry and the 

ferromagnetic quantum critical point 
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We study one-dimensional itinerant ferromagnets with Heisenberg symmetry near a ferromagnetic 
quantum critical point. It is shown that the Berry phase term arises in the effective action of itinerant 
ferromagnets when the full SU(2) symmetry is present. We explicitly demonstrate that dynamical 
critical exponent of the theory with the Berry term is z = 2 + O(e^) in the sense of e expansion, 
as previously discovered in the Ising limit. It appears, however, that the universality class at the 
interacting fixed point is not the same. We point out that even though the critical theory in the 
Ising limit can be obtained by the standard Hertz-Millis approach, the Heisenberg limit is expected 
to be different. We also calculate the exact electron Green functions G{x,t = 0) and G{x — 0,t) 
near the transition in a range of temperature, which can be used for experimental signatures of the 
associated critical points. 

PACS numbers; 71.10.Pm, 71.10.Hf, 75.10.Lp 



I. INTRODUCTION 

Quantum phase transitions from a paramagnetic to a 
ferromagnetic (FM) state in two and three dimensional 
itinerant electron system have been studied extensively in 
the pasfeii2i. The standard approach in studying such sys- 
tems, first advocated by Hertai, involves two steps. First, 
one decouples the electron-electron interaction term in 
the Hamiltonian using a Hubbard-Stratonovitch transfor- 
mation. Next, the fermions are integrated out to obtain 
an effective Landau-Ginzburg (LG) theory in terms of 
the FM order parameter which provides an effective de- 
scription of the transition. However, it has been recently 
pointed out that the procedure of integrating out the 
gapless fermions in the above mentioned scheme might 
be tricky^. In particular, it might lead to singular terms 
in the effective LG functional which can change the crit- 
ical property of the transition. 

In contrast, very little attention has been paid to such 
transitions in one dimension (ID), partly due to the Lieb- 
Mattis theoremi which states that the ground states of 
one-dimensional Hubbard models with nearest-neighbor 
hopping and density dependent short-range interaction 
are always spin singlets. However, numerical work car- 
ried out in Ref. Isl has shown that the presence of a next- 
nearest-neighbor hopping term in the above mentioned 
class of models can result in a FM ground state. The 
properties of such 'FM Luttinger liquids' has been stud- 
ied in detail in Ref. 0. Recently, a bosonization analysis 
has also been carried out to study the paramagnetic-FM 
transition in ID Hubbard models with Ising symmetry^. 
It is shown that in d = 1, the critical theory is below its 
upper critical dimension and is governed by an interact- 
ing fixed point. A similar conclusion was reached in Ref. 
1^ from analysis of rotor models of ferromagnets. 

In this work, we carry out a similar analysis which ap- 
plies for Heisenberg ferromagnets by explicitly keeping 
track of the full SU(2) symmetry of the model. In par- 
ticular, we show that at the Gaussian level, the presence 



of a Berry phase term in the effective action of such fer- 
romagnets does not change the dynamical critical expo- 
nent z of the transition: z ~ 2 + O(e^) in the e expansion 
sense. We point out that in spite of the same value of 
the dynamical critical exponent z at the Gaussian level, 
the universality class of this transition can be different 
from its Ising counterpart. Our analysis also shows that 
contrary to the claim of Ref. 0, the critical theory for 
Ising symmetry can be obtained by the standard Hertz- 
Millis procedure and in f D such a procedure is identical 
to the bosonization approach for Ising FM quantum crit- 
ical points. However, for the SU(2) case, it appears that 
the critical theory may not be described by the standard 
Hertz-Millis procedure. We also consider the effect of the 
FM critical fluctuations on the electron Green function at 
finite temperature. We obtain an exact expression for the 
electron Green functions G{x — 0,t) and G{x,t — 0) in 
a temperature range, which can be accessed experimen- 
tally by measuring the tunneling density of states (DOS) 
and the momentum distribution function respectively. 

The organization of the paper is as follows. In the next 
section, we obtain the effective action for ID quantum 
ferromagnets with Heisenberg symmetry at the Gaussian 
level and determine z by analyzing this action. Then 
we discuss the interacting fixed point in the e-expansion 
sense and the issue of universality of the Ising and Heisen- 
berg cases in the light of possible quartic couplings. This 
is followed by Sec. IIIII where we obtain the electron 
Green's function at finite temperature. We summarize 
and discuss our results in Sec. II VI 



II. CRITICAL THEORY IN THE HEISENBERG 
LIMIT 

We first derive the effective action for the critical the- 
ory by incorporating the transverse spin fluctuations in 
Sec. Ill AI This is followed by the analysis of the effective 
action to derive the dynamical critical exponent z in Sec. 
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Derivation of the effective action 



We start from a one dimensional t-J model given by 



n 



E 

<»i> 



(1) 



where ipi^j represents an electron at site i with spin cr, 
and S denotes the electron spin operators Sia = ipjo-a'ipi 
where Ca are the usual Pauli matrices. The correspond- 
ing action can be written as 



dx 



dT[lPHx,T)idr- fi)^ix,T)+H] (2) 



where r is the imaginary time and /? — l/ksT de- 
notes inverse temperature. Note that — and 

ip = {ipli^l)'^ spinor representations of the electron 
operators. We have switched to a continuous represen- 
tation of space and this amounts to replacing the last 
term in Eq. H by — / dxdx' J{x — x')S[x) ■ S(a;') where 
J((7) = Jcos(ga;a), where a is the lattice spacing and 
we shall consider J > throughout as appropriate for 
studying a ferromagnetic instability. 

Next, following Schultz^, we introduce new fermionic 
fields ip' defined as i}}{x^t) — U {x , T)ip' {x , t) . Here U 
represent SU(2) rotation matrices such that U'^{x,t)<t ■ 
n{x,T)U{x,T) = Gz- The unitary matrices U therefore 
correspond to the rotations of the local spin quantization 
axis of the electrons to z. The action S*, in terms of these 
new fields, becomes 

S = Sa + Sj 

So ~ dxdr -0'^ dr — fJ, — Cq — 2i cos {~idx — Cx) 4'' 



(3) 



where the fields C^, = -iU^di_,U = Y,a=x,y.z ^oS^y^ are 
the SU(2) gauge fields which describes the spin fiuctu- 
ations of the system, /i = (r, x) refers to the time and 
space components of the fields and 9^ — {idr,dx)- The 
gauge freedom here consists of rotation of the spin of the 
ip' fields about the z axis: tp' exp (i(T2<f>). Such a 

rotation changes fi^ — > 51^ + 9^$ and therefore leaves 
the action invariant. 

From now on we shall replace tp' by V' for notational 
convenience. The fermion fields can be written in terms 
of the right and the left movers: ipix) = i>L{x)e~'^^^^ + 
'!/'fl(a;)e''^^^, where kp is the Fermi wavevector. The right 
and the left moving fermions have energy dispersion 



Sj = - drdxdx' J{x-x')S'ix)S''{x') 



sgn{a)vFkx 



2m 



(4) 



where vp = 2ta sin {kpo) is the Fermi velocity and 
TO = 1/ [4ia^ cos (fci?a)] is the mass for the fermions and 
sgn(Q!) = +(— ) for a = R{L). In most of our analysis, 
we shall restrict ourselves to the linearized dispersion re- 
lation for the fermions, which amounts to neglecting the 
quadratic term in ea{kx). However, for deriving an ef- 
fective action for the spin-fluctuations we shall need to 
retain the quadratic term, and hence it can not be ne- 
glected at the outseliiS. 

In terms of these fields, the action reads 



S 

Sq 

Si 



S'o + 5*1 + 5*2 + 5*,/ 

d'kY,€ik)Go'4o.{k) 



(5) 



J d'kd^pJ2€ik+P) 

1pa{k) 



dk^ '"^^ 



s. 



Jd^k d^pY,4iik+p) 

4a{k) 



dkl 



Sj = - j d\Y.J^o.-{q)Sl{q)Sl,{-q) 



(6) 

(7) 
(8) 



notation for 
+ denotes anti- 



Here / d^fc is the shorthand 
P-^T.^J dkx/{2^), k^ = {iu;n,kx), 
commutator, a denotes the R, L indices for the right (R) 
and left(L) moving fermions, Gq^ = iujn — ea{kx) is the 
free fermion propagator, and Jaa'W) = J aa' cos.{qx) ■ 
We have neglected all 2k p terms in S since we shall 
be primarily interested in studying the action near a 
ferromagnetic instability. The matrix Jaa' — J + SJSaa' 
has an additional SJ term in its diagonal components, 
which is added to make the matrix J invertible and will 
be set to zero at the end of the calculation. The term S2 
Q represents the diamagnetic term of the SU(2) gauge 
fields and is known to be important for deriving a 
low energy effective description of the spin fluctuations. 
Note that to obtain this term, one needs to go retain the 
quadratic term in ea{kx)^- 

We then carry out a Hubbard-Stratonovitch transfor- 
mation to decouple Sj: 



S = 
Si = 

Srh — 



I d^kY,i,p^{k)si{-k) (9) 

-\ J d^kY,{J{q)VMq)K'{-q) (10) 



Notice that the field does not correspond to the phys- 
ical spin-density fields, but their conjugate. However, 
for studying the ferromagnetic instability, it is prefer- 
able to look at the physical spin-density fields. We 
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therefore introduce these fields p via a second Hubbard- 
Stratonovitch transformation which decouples 5*0: 

S = 5o + 5*1 + 5*2 + ^int + '5'hS 

Sint = I <fkY,iMk)[Sl,{-k) - p^(-k)] (11) 
5'hS = - J d^k'^Jaa>{q)pa{q)pa'{-q) (12) 

The next task is to integrate out the fcrmionic fields and 
obtain an effective low-energy description of the system 
in terms of the fields , p and 0. This procedure turns 
out to be quite subtle as discussed in Ref. 0, but can be 
carried ou^. The resulting effective action, to Gaussian 
order, reads 



S'nlsm 

^Bcrry 

S[p, 0] 



crry 

iV(0) 



(13) 
(14) 

(15) 



d\%{q)Y^p^{-q) 

a 

j '^^lYl i^'t^^^al'^t^a' - i(t)a{q)Pa'i-q)Saa' 

-Jaa'{q)pa{q)pa'{-q) (16) 



where we have expressed iSnlsm in terms of the unit vec- 
tor field n using the identitjiia: Tr [C^,f - Tr [a^C^ = 
{d^nf /4. The Berry phase term ^Berry , as is well known, 
can not be expressed in terms of n vector since fig de- 
pends on gauge choice of the SU(2) fieldsi^. The polar- 
ization tensor T\aa' in S[p, (p] is given by 



n-l,{q) = 5^o.'N{Q)- 



sgi\{a)vFq 



-iqQn 



sgn{a)vFq 



(17) 



Within the linearized dispersion and in one dimension 
Eq. ^] is known to be exacti^. 

We now introduce the spin-density fields in the sym- 
metric(S) and the antisymmetric(A) channels ps(A) = 
{pr + {~)pl), integrate out the auxiliary fields (pa from 
Eq. El Since the ferromagnetic instability corresponds 
to the instability of the spin-density fields in the sym- 
metric channel, we also integrate out the pA to obtain an 
effective action for ps 



ScS — S'nLSM + •S'Bcrry + S[ps] 

S'Bciry = J d'^q%{q)psi-q) 
S[Ps] = / d'^qps{q)Ssiq)psi-q) 



(18) 
(19) 



d^qps{q) 



^ r-hq']ps{-q) 



where r = (1 - JiV(O)), b = JN{Q)a?/2, and a is the 
lattice spacing. Here we have Wick-rotated back to real 
frequency and expanded J{q) ~ J(l — q^a^ 12) . Notice 
that the ferromagnetic instability in Eq.j^Olis signaled by 
r = or JiV(O) — 1 which agrees with the usual Stoner 
criteria for ferromagnets. 

We would like to stress that, unlike the Ising case stud- 
ied in Ref. 7, our final effective action fEa. ll8|) is not only 
S\ps\-, but also contains iSnlsm and ^Berry This is cru- 
cial for a correct description of the transverse spin fluctu- 
ations which are natural consequence of the SU(2) sym- 
metry of the problem. This point is best illustrated by 
studying the action (Eq. I18|l inside the FM phase. Here 
the spin-density fields are gapped out and (ps) = Af, 
where M is the magnetization of the system. In the FM 
phase, the low energy theory is therefore described by 
S'nLSM and SBcrry Notc that the presence of the Berry 
phase term is crucial for the fc? dispersion of the spin 
waves, as pointed out in Ref. [ij in the context of 2D 
ferromagnets. A similar analysis in this line for the Ising 
limit would yield Eq.l^Olas the final effective action which 
is exactly the same result obtained in Ref. X 

We also note that a similar procedure can be carried 
out in higher dimension. In that case, due to the pres- 
ence of a continuous isotropic Fermi surface, instead of 
two Fermi points, the index a becomes the momentum 
fc|| along the Fermi surface and the sum over a has to 
be replaced by an integral over k\^ . This changes the 

9o/(^f9^) term in 5s(g) (Eq. EOjl to |go|/(w/9) and im- 
mediately leads to z = 3. This result can also be al- 
ternatively obtained by finite dimensional bosonization 
techniquei^. 



B. Analysis of the effective action and critical 
theory 

We now analyze the critical theory for ID. Our starting 
point is Eg. 1181 derived in the previous subsection. From 
Eq. El the fluctuations of the n field are gapless in the 
Heiscnberg limit, it is not a priori clear whether the cou- 
pling of Ps to these fluctuations changes the universality 
class of the transition. To check this, we aim to integrate 
out the transverse fluctuation modes to obtain the flnal 
effective action in terms of the ps flelds. To this end, we 
rewrite the action S'nLSM and SBcrry using the CP(1) rep- 
resentation. In this representation, the unit vector field 
n is represented by two complex fields (zi , z ?) w ith the 
constraint \z\^ -I- |z2p = 1- Following R.efs. llMlfl . one 
can write 



Si 



NLSM 



dt dx ^ \ {df, - iAf_ 

.7 = 1,2 



)z,P + A(|z,f -1) 



where the U(l) gauge fields are given by 



(20) 



(21) 



(22) 
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The field A is a Lagrange multiplier field used to imple- 
ment the constraint of unit norm for the CP(1) fields. 
Note that in this representation, ^Bcny can be written as 



Sv 



(PqAo{q)psi-q) 



(23) 



and has the convenient interpretation of the action de- 
scribing coupling of a U(l) scalar potential to matter 
field. 

To describe the Berry term in a gauge invariant fash- 
ion, we now define a spin-current which satisfies the 
continuity relation 



dtPs + dxjl = d^J^{x) = 
and write the Berry term as 



BciTV 



(24) 



(25) 



Although this procedure seems ad-hoc, the additional 
term in S'serry introduced here could easily be obtained 
rigorously from our starting action (Eq. ^ if we carried 
out the analysis in the presence of a spin-current inter- 
action term i?currcnt ~ ^ J dx J dx' K{x — x')j(x) • j(a;') 
where j" = the spin-current operator. 

Such a current-current interaction term, although always 
present in principle, is usually neglected since it is very 
small {0{vp/c^)) compared to Sj. The only relevant in- 
formation, apart from the continuity condition 124f) . that 
we shall need for our analysis here is that the spin-current 
correlator in the symmetric channel is given hy^ 

^s{q) = -vl(^-l + Kcosiqa)NiO) 

(26) 

where r' — 1 — KN{Q). We therefore simply use this 
information and avoid carrying out a detailed analysis 
involving iJcuncnt- Note that the propagator ICs is always 
massive near the critical point. Hence it is not necessary 
to retain the q^ term in its expression. 

At this stage, we find that to understand the effect of 
S'scrry On the ferromagnetic quantum critical point, we 
need to obtain an effective action for the gauge fields 
Afj_. This involves integrating out the fields z and A from 
-Snlsm- It is well known that such an analysis can be 
reliably carried out for CP(iV — 1) theories only in the 
large N limi t^^'^^ . However, it is conjectured in Ref. 
that the qualitative results hold even for N = 2. W e 
therefore follow the large N analysis of Refs. Il5ll6l to 
obtain 



ScH — Sa-ausc 5'Bcrry 

+ s, 



y \ '-'matter 



(27) 



(fpA^{p) -p^iPu) Ay{-p) 



d^pA^{p)GlT{p)M-p) 



(28) 



where p^ = (po, vfPx)- Notice that the propagator G^^'^'^ 
is not invertible. This problem can be easily taken care 
of by simply adding a gauge fixing term, as is customary 
for the photon propagator in QED. However, we shall 
not need to invert the bare propagator at any stage of 
our analysis, and hence retain the form of G^^'"^ as given 
in Eg. 1281 The action S'mattor has now to be supplemented 
by a quadratic term in the spin current j| which involves 



S'mattcr = / d"^ q ps {q)Ss (q) PS {~q) 



+fAq)tCs{q)fA-q) 



(29) 



The final step is to integrate out the gauge fields. To 
do this, we first integrate out the matter field and 
obtain an effective dressed propagator cdrossod ^^^q 
gauge fields 



/ „2 1.2 I 1 
^dressed/j^\ _ _ / ^F^x ^ 45s (fc) 



LOVpkx 

1 

ilCsik) 



(30) 



We then replace G^''''' by G^''=''"'='^ in Eq.EHl and integrate 
out the gauge fields to obtain the contribution 6S of the 
gauge fields to the effective action for the matter fields 

T 18 

JlI — 



'if— 



5S 



5P{k) = 



•S'mattor -'r 5S 
j2 



d'J^{p){5P)^Ap)M-p) 

ICs{k)Ss{k) 
JCs{k)Lo^ + vlklSs{kf + l/^ 



(31) 



iK{k) 



-LOVpkx 



^\ -^vpkj vlkl + j^j 

Note that the term 6S couples the ps and j| fields. To 
obtain the effective action for the ps fields, we integrate 
out the jx fields and evaluate the resulting propagator 
at the critical point r = and retain the lowest order 
terms in lu and k. After some straightforward algebra, 
one obtains the correction term SS'[ps] 



SS'[pc 



d^kps(k) 



,2 1,2 



Ps(-fc)(33) 



where the ellipsis represent terms which are higher or- 
der in ijj and kx- So we conclude that the effect of the 
Berry term is merely to renormalize the coefficients of the 
existing terms in S[ps\- The critical theory obtained by 
integrating out the transverse fiuctuation at the Gaussian 
level thus has z = 2. By the usual e expansion argument, 
we expect this analysis to give correct value of z to O(e^). 

Next we comment on the universality class of the 
transition. To determine the universality class, since 
we are below the upper critical dimension of the the- 
ory, we need to retain the possible fourth order dia- 
grams while integrating out fermions in Eq. 1111 Two 
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a 

z 


z 


A" 

a 


V 



a. 



n 



quires a consideration of both coupled longitudinal and 
transverse modes at the same footing. In this work, we 
shall be content with pointing out the necessity of such 
an analysis. 

The quartic terms mentioned above makes the Gaus- 
sian fixed point, described by Eg. 1181 unstable. However, 
all such diagrams, including the Ising type quartic terms, 
are non-zero only when the fermionic dispersion has a fi- 
nite curvatur e ''i^^ . Thus, above a temperature scale set 
by this curvature, the properties of the system is expected 
to be well described by the Gaussian theory described by 
Eq. El We discuss this point in more details in the next 
section. 



III. ELECTRON GREEN'S FUNCTION 



FIG. 1: Representative diagrams for quartic coupling terms 
for Heisenberg theory. The left diagram is also present for 
Ising ferromagnets and generates the usual pg interaction 
term. The diagram on the right side generates quartic cou- 
plings between ps and n fields (see discussion in text), and 
has no analog in Ising ferromagnets. The internal lines in- 
volve Fermion Green functions and the vertices represent ei- 
ther az or A^^ — aadG^^ /dkfj, as shown. Both the diagrams 
are nonzero only if the Fermion dispersion has a non-zero cur- 
vature. 



such representative diagrams are shown in Fig. ^ The 
diagram in the left panel is present for the Ising fer- 
romagnets studied in Ref. 0, while that in the right 
panel is a consequence of the Heisenberg symmetry of 
the problem. After some algebraic, one can show that 
the diagram in the left panel of Fig. ^ generates Ising 
type contributions ~ J dxdtpg while the diagram in 
the right panel generates, for example, quartic couplings 
~ J dxdtp^{dx'n.)^ (apart from other similar quartic and 
higher order terms). The former type of terms are exactly 
those that are expected at quartic order for the Ising fer- 
romagnets, either from a bosonization analysis^ or Hertz- 
Millis approacbi*^. Therefore, our analysis shows that 
the standard Hertz-Millis procedure in ID, gives exactly 
the same critical theory as obtained in Ref. 7 by bosoniza- 
tion for the Ising case. 

The latter type of diagrams are absent in the Ising 
model studied in Ref. 0. In the FM phases, where the 
amplitude mode is gapped, these terms represents a triv- 
ial renormalization of the parameters such as spin-wave 
stiffness of ^nlsm (Eq- IT^ . Near the transition, they 
generate quartic coupling between the amplitude and the 
spin-wave modes. Notice that these terms are different 
from the standard m'* quartic terms expected from the 
usual Hertz-Millis analysis of a Heisenberg magnet with 
a vector order parameter. We expect such quartic cou- 
pling terms to change the universality class of the transi- 
tion from that of Ising ferromagnets. A determination of 
the universality class would therefore require a detailed 
analysis of Eq. E| supplemented with all such possible 
relevant quartic coupling terms. In particular, this re- 



It is well known that the presence of a quantum criti- 
cal point shows up in experimentally accessible quantities 
such as tunneling DOS and the momentum distribution 
function of the electrons. To obtain these quantities one 
needs to estimate the effect of critical fluctuations on the 
electron Green function. However, the fixed point con- 
cerned here is an interacting one. This means to compute 
the electron Green function at T = near the critical 
point, one needs to solve the problem of electrons cou- 
pled to interacting bosons. For the Ising ferromagnets 
these bosons can be described by an action 



S 



dxdr 



(_a2 _ ^g2 ^ f^g4^ ^ ^ m(9,0)^ (34) 



where u is determined by the curvature of the fermionic 
dispersion. Obtaining the electron Green function cou- 
pled with these interacting bosons, even for the Ising 
ferromagnets, is a difficult task. However, it might be 
still be possible to compute the properties of Green func- 
tion at finite temperature using the Gaussian fixed point, 
if the RG flow, which can be computed for Ising ferro- 
magnets from Eq. 1341 generated by temperature takes 
us to a regime where bq^ > r{T) 3> u*, where qt = 

y^T/Tp kp / [1 — r(T)]^^^ is the thermal wavevector and 
u* is the value of the coupling u at the interacting fixed 
point. 

To explain the last sentence a little better, we consider 
the RG fiow of the action S (Eq. I34|l at finite temper- 
ature, to one loop. As noted in Ref. 7, the flow of u 
and b are negligible as r flows away from the fixed point 
according to 



dr 
'dl 



2r 



3" I ^2 
TT 



- coth — 
2/ \2T 



(35) 



where T = Toexp(zZ) is the scaled temperature. There- 
fore we may envisage that when e is small, there will be 
a temperature Ti, above which r(T) m*, so that the 
effect of interaction can be neglected. The key question is 
then whether at and above Ti, r(T) is small compared to 
the term in S (Eq. I34|l . A simple estimate shows that 
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for this to happen, we need to be below a temperature 
T2 given by 



T2_ 

Tf 



r{T2 



(36) 



For T > T2, r can not be neglected and the behavior of 
the system is similar to that of a Luttinger liquid at finite 
temperature. 

If indeed Ti < T2, there exists a finite window where 
the effect of the critical fluctuations on the electron Green 
function can be calculated neglecting the effect of inter- 
action. We refrain from estimating Ti and T2 since they 
need the knowledge of precise value of u* and is therefore 
non-universal in the RG sense. Instead, in the rest of this 
section, we shall assume that such a window exists and 
compute the behavior of the electron Green function for 
Ti < T < T2. Note that such an analysis, which requires 
only properties of the Gaussian fixed point, also holds for 
the Heisenberg ferromagnets, since the Gaussian action 
is the same in both cases. The question of existence of 
such a window of course has to be separately investigated 
for the Heisenberg ferromagnets. Such an investigation 
requires a detailed analysis of quartic couplings for the 
Heisenberg theory and is beyond the scope of the present 
work. However since all such quartic terms depends on 
the curvature of the Fermionic dispersion, we expect such 
a window to exist also for the Heisenberg ferromagnets 
also when the curvature is small. 

To calculate the Green function for the electron we 
start from the action 

S = Sq -\- Si 

50 = f d^Qd^k^Uk + q) 

51 - ^d'#a(Q)FRiiA('/)0a(-g) (37) 

a 

where (j) denotes the bosonic fields with the propagator 
^RPA given by 



^rpa(<7) 



Jiqx 



l + ^fe)Eana^(9) 



(38) 



where H^^^ — N{T)avFqx/ {i^n — ctVFQx) is the finite 
temperature generalization of H~^, (q) in Eq. 1171 for 
q,/kF < 1, and N{T) = Y.kdf{^)/d{e) ~ iV(0) is the 
free fermion DOS. Note that the starting action given by 
Eq. 133 can be easily obtained from Eos. II II and 1121 with 
the propagator for the p{q) fields replaced by Frpa— . 

To obtain the Green's function of the electrons, we 
proceed following the analysis carried out in Ref.^3 The 
equation for the Green function is given by 

S{k-k/) = Y.['^o\k)Sik-ki)S,,, 

ki 

+iia,)aa'Mk - fci)] C'/ik, k') (39) 



which in ID can be solved exactly using Schwinger 
ansataiSiSS. After some straightforward manipulations, 
we get, G"'^ {x,t) = Scra'G{x,T), where 



Gix,T) 

Ga{x,T) 

Goa(a;,r) 



^e'"'=--G„(x,T) 

a 

Goa{x, t) exp [Q{ax, t) 



cPk g^C^^^-'^"'^) 

iLO„ 

-1 



(40) 



avpk^ 



27r J sinh [(ax — ivft) / ^] 



(41) 



Q{ax,T) = cfqFRPAiq) 



1 — cos (qxX — f2„r) 

(ifl„ - avFqxf 

(42) 

Here (£!„) denote fermionic (bosonic) frequencies sat- 
isfying the usual anti-periodic (periodic) boundary con- 
ditions, ^ = hyF/ksT is the thermal correlation length, 
and one needs to make the standard Wick rotation t = it 
to obtain the Green functions in real time. 

To proceed further, we note that experimentally ac- 
cessible quantities such as the momentum distribution 
7ia{k) or the tunneling DOS p{lo)^ do not probe the full 
Green function Gq(x, t), but only Gq(0, t) and Ga{x, 0), 
since 

/dk 
exp{ikxx)G aix,0), (43) 

p{uj) = / dtIm[exp(-iwt)G„(0,i)] . (44) 



Therefore we resort to the simpler task of computing 
G(x,0) and G(0, r). The frequency sums in the expres- 
sions of Q{x,0) and Q(0, r) (Eq. 021) can now be evalu- 
ated in a straightforward manner. We express the final 
result in terms of the crossover temperature Ti and the 
corresponding length scale ^1 — hvF /ksTi = k^^TF /Ti. 



Q{x,Q) 



exp 



, 1 — cos(g'x') 



Q(0,i) 




exp 



X [1 + 7(17')] coth 



^-'^°^('^'^')cothf|- 



where all momenta, length, and time scales are ex- 
pressed in terms of q' = qx£,i, x' = x/^i, s = T/Ti is 
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FIG. 2: Plot of G{x, 0) (dashed line) and G(0, t) (solid line) 
near the critical point (r = 0) at T = Ti . The decay of G(0, t) 
is slower than G{x, 0) as discussed in the text. 



the dimensionless temperature, and t' = vpt/S,!- Here 
/i((7') = vp{q')/vF is the ratio of the renormalized and 
the bare Fermi velocities and 7(9') is a dimensionless pa- 
rameter given by 



iW) = 



1 + /i(g') 



l\2 



- 1 



2/i(g') 



(47) 



and = qc^i is the dimensionless infrared cutoff. Note 
that in our notation, the non-interacting limit corre- 
sponds to ^{q) = 1 = r. 

Before looking at the Green function in its full gen- 
erality, it is instructive to check that it reproduces the 
well- known result for the Luttinger liquid Green func- 
tion when &' = and r 7^ 0. In this limit the integrals 
in Eq. 0^1 and 0^1 can be calculated analytically and one 
gets at T = Ti, 



G^^(x',0) - 



Gj;^(0,t') = 



-i\ ( M(0)/ei 



27r / \sinh(Q!a:/^i) 

^sinh(|xV(0)/Ci) 
Ga{ax' = i',0). 



7(0) 



(48) 



From Eq. 0H1 one can identify 7(0) to be the anoma- 
lous dimension^^. Note that, as is well known for Lut- 
tinger liquids, the Green functions G{x,Q) and G(0,t) 
has identical behavior. This is a consequence of the lin- 
ear dispersion of the bosonic density fluctuations. As 
we approach the critical point, r — > 0, and consequently 
one must retain the h' term, in the expression of 7(17') • 



When 6' 7^ 0, it is not possible to evaluate Q analyti- 
cally at finite temperature. However, since the electrons 
now see the critical bosonic fluctuations with dispersions 
Lo ~ q^, we expect the Green functions G+(a;,0) and 
G(0,i) to have different behavior. A plot of the Green 
functions shown in Figs.|21 conflrms the above qualitative 
discussion. The characteristic scale of decay for G{x, 0) 
is seen to be much shorter compared to G(0, t). This fea- 
ture arises from the presence of critical fluctuations with 
uj ^ q^ dispersion ( in contrast to w ~ g^, fluctuations in 
standard Luttinger liquid ) and is therefore a signature 
of the quantum critical point at finite temperature. We 
expect that a measurement of p{oj) and n{kx) will probe 
this behavior. 



IV. DISCUSSION 

In this work, we have shown, by explicitly keeping 
track of the full SU(2) symmetry of the problem, that 
in ID the Luttinger liquid to FM transition for Heisen- 
berg ferromagnets has z = 2 -I- O(e^) in an e expan- 
sion sense. The analysis done here was carried out at 
a Gaussian level, which is exact in ID as long as the 
dispersion of the Fermions are linear. The curvature in 
Fermionic dispersion introduces quartic coupling terms 
beyond the Gaussian approximation. Although we have 
not carried out a full analysis of the problem with such 
coupling terms in the presence of the SU(2) symmetry, 
we have identified the typical quartic terms and discuss 
their consequences. In particular, we pointed out that for 
a ferromagnet with Ising symmetry, our analysis reduces 
to the standard Hertz-Millis theory. The Hertz-Millis ap- 
proach is therefore coiirpletely equivalent to the bosoniza- 
tion approach of Ref. 7 in ID for Ising ferromagnets. In 
contrast, we find that for the Heisenberg ferromagnets, 
our theory does not reduce to that obtained by the tra- 
ditional Hertz-Millis approach. We identify that this is 
due to additional coupling terms between the longitudi- 
nal and the transverse modes, shown in right panel of 
Fig. n which has no analog in the Ising case and are also 
different from the usual expected m"* interaction terms 
that one obtains using the Hertz-Millis approach. It will 
be interesting to carry out our analysis in higher dimen- 
sion and compare the results with those of Ref. 0, where 
a possible problem with such quartic terms in the 
usual Hertz-Millis approach has been discussed from a 
different point of view. A detailed analysis of this issue 
is left as a subject of future work. 

We have also obtained an exact expression for the elec- 
tron Green functions G{x^ 0) and G(0, t) for a range of 
finite temperature. Although we have not obtained the 
expression of the Green functions at the interacting fixed 
point, we hope that over a range of temperature where 
the properties of the system is expected to be well de- 
scribed by the Gaussian fixed point, G{x,Q) and G(0,t) 
can be probed by experimentally accessible quantities 
such as the momentum distribution and tunneling DOS 
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of the electrons. 
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